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^ Abstract 

^ We use a method recently suggested for evaluating the slope of the Isgur-Wise 

function, at the zero- recoil point, on the lattice. The computations are performed 
in the quenched approximation to lattice QCD, on a 24^ X 48 lattice at /3 = 6.2, 
using an 0(a)-improved action for the fermions. We have found unexpectedly large 
finite-volume effects in such a calculation. These volume corrections turned out to 
be purely geometrical and independent of the dynamics of the system. After the 
study of these effects on a smaller volume and for different quark masses, we give 
approximate expressions that account for them. Using these approximations we find 
^'(1) = —1.7 2 and ^'(1) = —1.4 _^ for the slope of the Isgur-Wise function, for two 
mesons composed of a heavy quark slightly heavier and lighter, respectively, than the 
charm quark, and in both light antiquark whose mass is about that of the 

strange quark. 
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1 Introduction 



In recent years there has been an increased theoretical and experimental interest in the 
study of decays of hadrons containing a b quark. Due to the appearance of a spin-flavour 
symmetry in the dynamics of heavy quarks ([1]), systems involving a heavy quark are in some 
ways simpler than those involving only light quarks. The study of these systems has several 
objectives, including the extraction of the elements Vet and Vub of the Cabibbo-Kobayaslii- 
Maskawa matrix, a determination of the applicability of the "Heavy Quark Effective Theory" 
(HQET) predictions for the physics of b and c quarks, and derivations of bounds for the effects 
of physics beyond the standard model. 

Recently the ARGUS [2] and CLEO [3] collaborations have measured the semileptonic decay 
B D*W, and have presented results for ll/cfcl^l'^), where ^{lu) is the Isgur-Wise function, 
C{uj), up to short-distance and heavy-quark mass corrections [4] and ui is the dot product 
of the four-velocities of the B and D* mesons. The Isgur-Wise function is the unique form- 
factor for all semi-leptonic B — ^ D-decays in the limit of infinite b and c quark masses. The 
extraction of the \ Vcb\ matrix element requires a precise knowledge of the Isgur-Wise function 
in a region of values of to close to the zero-recoil point (i.e. lu = 1; the Isgur-Wise function 
is normalised such that ^{1) = 1). 

Lattice computations allow for the determination of the Isgur-Wise function for discrete 
values of (J > 1 ([5] - [6]) and these can be used to reconstruct the Isgur-Wise function near 
to the zero-recoil point, and in particular to estimate its derivatives at lu = 1. Although 
much of the systematics is understood in this kind of calculation, [7], there are still some 
problems which lead to uncertainties in the determination of ^{lu) near to the zero-recoil 
point (extrapolation of the results from ui > 1 to ui = 1, determination of the possible ui- 
dependence of the systematic effects which may distort this extrapolation, etc..) and could 
lead to uncertainties in the extraction of 1 • For this reason, the authors of [8] suggested 
a new method, in which the slope and higher derivatives of the Isgur-Wise function at 
uj = 1 are computed directly on the lattice, hence avoiding any kind of extrapolation (which 
necessarily relies on a certain parametrization of the (j-dependence of Isgur-Wise function) 
from a; > 1 to the zero-recoil point. The results obtained by using this new method would 
complement those obtained by using the traditional one ([5] - [6]) and would help to unravel 
possible sources of systematic errors in lattice calculations. 

In future it is likely that it will also be possible to study the Isgur-Wise function near the 
zero-recoil point by simulations using the HQET. Results from an early study can be found 
in ref. [9]. 



1 



In [8] it was proposed to study correlators of the type: 

ci'\t) = -\j:xi{o\jp{x)4{om, (1) 

X 

cf(t,g = -\j:4{o\My)v\x)jUo)\o), (2) 

in addition to the usual two- and three-point correlators, {C2^\t) and C^^\t)) given by: 

Ci'\t) = X:(0|Jp(^)4(0)|0), (3) 

X 

Ci'\t,t,) = j:iO\Jpiy)V\x)jUom, (4) 

where Jp{Jp) is an interpolating operator which can annihilate (create) the pseudoscalar 
meson P, and represent the fourth component of the vector current Qj^Q. In [8] it was 
shown that, from the measurement of these correlators, the slope of the Isgur-Wise function 
at the zero-recoil point could be extracted. Actually, if local interpolating operators are used 
and in the infinite- volume limit, it can be shown that: 

1 '1 ■ ''^"(l)) = Rf>-Rf\ (6) 



2m2^2 dcu 

where ^{co'), as we mentioned above, is the Isgur-Wise function up to short-distance and 
heavy-quark mass corrections and the ratios ^ defined by: 

^2,3 = -(5y- W 

The second derivative of the Isgur-Wise function would be determined by a similar set of 
correlators to those of eqs.(l-2) in which has been replaced by x^, and in principle, higher 
order derivatives could be determined alike [8]. 

Below we present the results of an exploratory study to check the feasibility of the method 
proposed in [8]. The main result of this paper is that in the evaluation of correlation functions 
of the form in eqs.(l) and (2) there are large finite- volume effects. Although these effects, 
for a fixed meson mass and time, decrease exponentially with the length of the lattice in 
the spatial directions (with a rate given, to a good approximation, by [m/t]^)^ they are 
large on currently available lattices. The origin of these effects is purely geometrical and 
understandable, and they do not depend on the dynamics of the system. The ordinary 
two- and three-point correlators (without the insertion of the x\,x\, ... terms) used in lattice 
calculations are free of this new source of systematic errors. The problem arises because 
when the terms included in the definition of the correlators, the lattice theory 



2 



becomes sensitive to momentum contributions in addition to the unique momentum mode 
which determines the ordinary two and three point correlators^. The behaviour of these 
geometrical volume effects, and in particular, the dependence on the mass and time of 
evolution, is different from that in usual correlation functions. Similar effects Avill be present 
in any direct computation of a derivative of a matrix element with respect to an external 
momentum. Since the origin of these effects is purely kinematical, they can be controhed, and 
below we present a set of approximations which will allow us to extract the slope of the Isgur- 
Wise function at zero-recoil in simulations on finite volumes. Using these approximations 
we find ^'(1) = —1.7 ^^'^^ ^'(1) = —1.4 for two mesons composed of a heavy quark 
slightly heavier and lighter , respectively, than the charm quark, and in both cases, a light 
antiquark whose mass is about that of the strange quark. These values agree well with those 
obtained in the lattice computation of ref. [7]. In this reference, for the same heavy and 
light quark masses and using the traditional method of extrapolating from a; > 1 to the 
zero-recoil point, were found the following values for the slope of the Isgur-Wise function: 
^'(1) = -1.4 +J and ^'(1) = -1.4 respectively. 

The remainder of this paper is organized as follows. In Section 2 we will present and study 
the origin of these new volume effects in the simplest case: two-point function correlators 
with local interpolating operators Jp. The results of this section will permit us to study in 
Section 3 the case of extended interpolating operators and the computation of the slope of 
the Isgur-Wise function at zero-recoil. Finally, in Section 4 we will present our conclusions. 



2 Geometrical Finite- Volume Effects 

In this section we demonstrate the existence of finite-volume effects in correlation functions 
of the form in eq. (1), and show that these effects are large on currently available lattices. 
Consider the evaluation of such a correlation function on a N"^ x Nt lattice, of lengths 
T and L in the temporal and spatial directions respectively, and with lattice spacing a 
{L = Na,T = Nto): 

C?\t) = -\T.4{0\M^)4mO), (7) 

X 

where Pmin = 27r/L, Z^(p^) = |(0| Jp(0)|]5, ^4)1 and the sum X)a is over all the excited states 

^This momentum mode coincides with the external momentum provided in the definition of the 
correlators. 
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A that can contribute to the correlation function. rriA is the mass of state A^. In general 
Z(p^) depends on the momentum \p\ if extended interpolating operators Jp are used. 

In the infinite-volume limit L, T — ^ oo, and with t ^ I so that only the ground state 
contributes significantly^, the correlation function becomes [8]: 

ci'\t)^Cf\t)xi-^^-±- + 2R^^\ (9) 

where 

Ci'\t) = ^e-\ (10) 

We now study the finite- volume corrections to eq. (9). For purposes of illustration let us take 
the interpolating operators Jp{x) to be local ones, so that the corresponding wavefunction 
factor Z(p^), in eq. (8), is independent of momentum, and can be taken out of the sum over 

(2) 

n, and also i?y = 0. Having determined the mass of the ground state from the behaviour 
of the correlation function C2^\t) with time (see eq. (10) ), the ratio of C2'\t) and 6*2°^ (t) 
can be calculated exactly by computing the discrete sum in eq. (8): 



. ||, (13) 



1 " 2 1 2 p-(\/m2+p^.^n2-m) t 

2 E J: ^=__^e— . (14) 



Assuming that: 



i) The lattice spacing in the spatial directions is zero, but that L is finite, so that the 
sum over can be replaced by an integral, i.e. 



I 1 



which implies that the sum over n in eq. (14) must be extended between — oo to -|-oo. 



^ There may also be contributions from multiparticle states, but we will shortly restrict the discussion to 
the ground state only. 

''In the following, except where expHcitly stated, we will neglect the contributions from the excited states. 
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ii) L 3> 1 and therefore we can perform an asymptotic expansion in powers of 1/L. 



the result of eq. (14) can be approximated by: 



R2\t) ~ ( - At - — ) X I 1 - — f . I • (16) 

|^ + ife)^ sinh(I/(^ + ^)^^ ' 

In general eq. (16) is a good approximation to eq. (14) when both ■n/mL and ■n'^t/mL?' 
are small quantities; the higher the mass m, and the smaller the time t, the better is the 

('21 

approximation. In Figs.l and 2 we test the validity of eq. (16) by plotting i?2 {t) as a 
function of the time, t, for different volumes and masses of the scalar particle. The plots 
correspond to masses given by ma = 0.3 and 0.6, and volumes V = 24? and 36^. In each case 

(2,) 

the dashed line represents i?2 {t) computed numerically from the exact expression given in 
eq. (14), the sohd line corresponds to the infinite-volume result, R^\t) = (— l/2m^ — t/2m), 
and the squares stand for the approximation of eq. (16) to i?2 • The first observation which 
can be made from these plots is that there is a considerable difference between infinite- (solid 
hue) and finite- (dashed hue) volume predictions. A second observation is that eq. (16) is in 
most cases a good approximation to the finite- volume prediction of eq. (14), i.e. for a wide 
range of volumes, masses and times, the squares follow the dashed line. 




Figure 1: Comparison between different approximations to i?2 {t) (see text), a) m=0.6; b) 
m=0.3. In both cases V = 24p. 
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Figure 2: Comparison between different approximations to i?2 (t) (see text), a) m=0.6 ; b) 
m=0.3. In both cases V = 36^ . 

From eq. (16) we see that the finite-volume eff'ects decrease exponentially with L and, for 
sufficiently large times, they are given by: 



where the dependence on L, m and t is explicitly shown. 

The origin of these geometrical volume corrections is due to the fact that more than one 
mode contributes to the sum over n in eq. (8). By contrast, in the usual two-point functions, 
C2°^(t), only the term with n = contributes in the sum, and thus this type of geometrical 
finite- volume effect is not present. However, both C2\t) and C2'\t) suffer from dynamical 
volume effects. These effects are either due to self interactions of the particle which is being 
studied with its infinite number of mirror copies (in simulations with periodic boundary 
conditions) [10], or due to the modification of the hadronic wave function in situations for 
which the physical volume of the lattice is comparable to the size of the hadron [11]-[12]. 
For small lattices, the latter type of corrections is the dominant effect, whereas for physical 
volumes of the order of 1.5 fm the predominant effect is due to the interaction of the particle 
with its mirror copies ([10], [12]). However, all of these corrections are small for the masses 
and volumes considered in this paper, and in what follows they will be ignored. 

As was mentioned in the introduction, the slope of the Isgur-Wise function in the infinite- 




(17) 
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volume limit can be extracted from the correlators €'2 (t) and C3 ' (t). For finite volumes such 

(2') 

a determination contains similar corrections to the ones studied here for C2 (t), corrections 
which decrease exponentially with the length of the lattice. The slope of the Isgur-Wise 
function at a; = 1 can also be estimated by computing the finite difference ^{idmin) — 'C(l)) 
where cdmin is the discrete value of cu nearest to a; = 1 available on the lattice being used. 
^'(1) derived in such a way suffers from finite- volume corrections of the order of 0{ujrnin) 
or eciuivalently of the order of 0(1/L^). These finite- volume corrections, therefore, formally 
decrease more slowly Avith the volume of the lattice than those affecting C2'\t) and C^^\t). 

In the remainder of this section we present the results from numerical simulations using 
local interpolating operators. We have measured R2 (t) for a heavy meson with a mass of 
about 1.6 GeV in a study on a 24'^ x 48 lattice at 3 = 6.2 and using the 0(a)-iniproved 
quark action (SW or "clover") proposed by Sheikholeslami and Wohlert, ([14, 15])^. The 
meson was composed of a heavy quark with a mass corresponding to Kh = 0.133 and a 
Ught antiquark with ki = 0.14144 ^. The mass of the meson was found (from 60 gauge field 
configurations) to be ma = 0.598 [13]. 

In Figure 3 we compare the prediction of eq. (14) with the data. The solid line corresponds 
to the infinite-volume prediction R^\t) = (— l/2m^ — t/2m), the dashed line to R^\t) 
computed using eq. (14) ^ (including volume effects) and finally the squares stand for the 
measured data (from 10 gauge field configurations) on the lattice. From Figure 3 it is clear 
that the dashed fine fits the data much better than the solid line, demonstrating that the 
discrepancy between the infinite-volume prediction (9) and the data is largely due to the 
geometrical finite-volume effects. However, even the dashed line does not reproduce the 
data perfectly, and we understand the small discrepancy between the two as being due to 
discretisation errors (i.e. errors due to the finiteness of the lattice spacing)^. This will be 
discussed further below. We have checked that the small discrepancy between the dashed 
fine and the data is not due to the contribution of excited states, since the size of such a 
contribution would be incompatible with the behaviour of the two-point function C2'\t). 

We have also computed the ratio R2 (t) on 36 gauge configurations, for a light meson 
composed of a degenerate (light) quark - anti-quark pair, in a simulation on a 16^ x 32 

^For general details of the simulation, which was performed on the 64-node i860 Meiko Computing Surface 
at the University of Edinburgh, see references [13] and [16]. 

^In this simulation the hopping parameter of a massless quark is Kcrit = 0.14315(2), and that of the 
strange quark is = 0.1419(1) ([13]). 

'^We have appropriately modified eq. (14) in order to include the contribution of the meson propagating 
backwards in time. 

^The 0(a)-improved action which has been used in this work is free of niQU discretisation errors. The 
leading discretisation errors are of order amga or {mQof' ([15]) which for a heavy quark may be of a 
comparable size, mg is the mass of the heavy quark. 
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Figure 3: Local interpolating operators: Comparison between the approximation of eq. (14) 
(dashed hne), the infinite- volume prediction (sohd fine) and the measured data (squares). 
Kh = 0.133, Ki = 0.14144, (3 = 6.2 on 24^ x 48 lattice. 

lattice at = 6.0, again using the S'l/V'-improved action for the fermions^. The mass of the 
light quark corresponds to k = 0.1425 (and that of the corresponding meson Avas found to 
be ma = 0.448 ^^): in lattice units the quark mass is about four times smaller than the 
heavy quark mass in the simulation described above. Thus the discretisation errors should be 
reduced considerably Avhile the finite-volume effects are going to be even more pronounced, 
and indeed Ave now find that this is the case. In Figure 4 we compare the prediction of eq. (14) 
with the measured data. The meaning of the different curA'es is the same as in Figure 3. It 
can be seen from Figure 4 that the agreement between our theoretical prediction and the 

(2) 

data is excellent, in spite of the fact that the finite-volume corrections in the ratio i?2 (t) 
are very large for such light mesons. 

We now try to justify our hypothesis that the small discrepancy between the theoretical 
prediction (dashed fine) and data (squares) in Figure 3 is due to discretisation errors. We 
have tried to reduce these errors by replacing the continuum energy-momentum dispersion 
relation in eq. (14) by the lattice dispersion relation of a free boson and the residue factor 
1/2E by l/2sinh(£') ([18]). More precisely we have replaced in eq. (14): 

g-(V™HpL„"=-™) * ^ e-^^'^^^™''^™^""')-™) *, (18) 

^This simulation was performed on the CRAY-YMP at Cineca in Bologna, and details of the simulation 
can be found in ref.[17]. 
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Figure 4: Local interpolating operators: Comparison between the approximation of 
eq. (14) (dashed hne), the infinite- volume prediction (sohd fine) and the measured data 
(squares). Ki = 0.1425,/? = 6.0 on 16^ x 32 lattice. 

1 ^ ^i"h(ma) 

V^l + PLnnVrn" sinh(Ei^(m2,pLrx«')a) ' 

where 

E'^^im^q^) = - arcsinh (Jsinh^(— ) + ^in\^)] . (20) 

Having made these changes, we again compare our prediction, now given by eq. (14) and 
eqs. (18-20) , with the data for the ratio R2'\t) of Figure 3. The result is plotted in Figure 5. 
The meaning of the different curves is again the same as in Figure 3. The agreement of the 
predictions with the data is excellent. This fact can be further underlined by noting that if 
we fit i?2 {t) to the theoretical prediction with the mass left as a free parameter we obtain 
ma = 0.60(1), in good agreement with the value of 0.598 obtained from C2°^(t). The 
result of such a fit is ma = 0.63(1) when the substitutions of eqs. (18-20) are not taken into 
account. 

In the light-light case of Figure 4, the inclusion of finite-lattice-spacing effects leads to minor 
changes, and a slightly improved i^i the comparison of the ratio R^\t) with the theoretical 
prediction. 

In summary, we repeat that, by making a reasonable ansatz for the discretisation errors. 
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Figure 5: Local interpolating operators: Comparison between the approximation of eq. (14) 
and eqs. (18 - 20) (dashed hne), the infinite- volume and zero lattice spacing prediction (solid 
fine) and the measured data (squares). Kh = 0.133, ki = 0.14144, /3 = 6.2 on 24^ x 48 lattice. 

we can account for the small discrepancy between the theoretical finite-volume prediction 
and the data in Figure 3. Moreover, for light mesons (Figure 4), there is no corresponding 
discrepancy. These observations, and the fact that the discrepancy cannot be explained 
by contamination by excited states, leads us confidently to attribute it to the presence 
of discretisation errors. We stress, liOAvever, that the main point of this section was the 

(2) 

demonstration that the infinite-volume prediction for i?2 has large, controllable, finite- 
volume elfects, as is clear from the difference between the solid and dashed lines in figures 
1-5. 

3 The Slope of the Isgur-Wise function 

In this section we will extract the slope of the Isgur-Wise function by computing two- and 
three-point correlation functions, using extended (or "smeared") interpolating operators in 
order to enhance the contribution of the ground state. In this exploratory study we restrict 
the computation to two values of the mass of the heavy quark, and to a single value of the 
mass of the light quark (as explained below). 
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Following the ideas of [8], let us consider now^°, 

ci'\t,ty) = -lj:xi{o\jp{y)v\x)juom (2i: 



mm) 



^ ^inprnin^S (22) 

where represents the fourth component of the vector current ZyQ^^Q. Zy is the renor- 
malisation constant relating the lattice vector current to the physical one. is the fourth 
component of the matrix element defined by 

M^{n) = {P{p')\Qi^Q\P{p)) = {p + p'n{^). (23) 
where p' = 0, p = (0, 0, npmin) and 

e» = (l + /3(a;) + 7(a;))e(a;), (24) 

^ = l-^T^' (25) 
q = p-p'. (26) 

/3(li;) represents the radiative corrections and 7(^1;) corresponds to corrections proportional 
to inverse powers of the heavy-quark mass (/3(1) = 7(1) = for the degenerate case we are 
considering here [1], where by degenerate we mean that the initial and final state mesons 
have the same masses); in what follows we will neglect the corrections j{uj), which appear to 
be very small in the range of masses studied in this work [7]. ^{ui) is the Isgur-Wise function 
and Lo is the dot product of the four velocities of the initial and final state mesons. The 
Isgur-Wise function is normaUsed to 1 at the zero-recoil point, i.e. ^{1) = 1. In eq. (22) 
it has been assumed that t 3> 1 and tj, — t 3> 1 and that therefore only the ground state 
contributes to the correlation function. 

In the limit L — ^ 00 the correlation function behaves like^^: 



where 



Ci'\t,) = Zy?^e-^'^M\0), (28) 
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For the present we will ignore possible discretisation errors. We will come back to this point later. 



^^Note that the ratio R^^\t) does not depend on the renormalisation constant Z 
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M^'{0) _ 1 4 




(29) 



■M 



M4(0) 2m2^2 duj 




d^ 



d^ 




(30) 



du) 



dM{p^) 
dp"^ 



(31) 



and was defined in eq. (11). The idea of the authors of ref.[8] was to obtain the 
derivative of the renormalisation-group-invariant Isgur-Wise function ^{uj) at the zero-recoil 



As was pointed out in [8], this method allows for a direct determination of the slope of the 
Isgur-Wise function at zero-recoil, without the need for an extrapolation of the computed 
values of ^{lo) to o; = 1. Moreover the derivatives are computed in terms of ratios of 
correlation functions, and it is likely that there is some cancellation of systematic errors and 
statistical fluctuations in these ratios. 

We have measured the ratios R2 3 using 49 gauge field configurations in our simulation on 
the 24'^ X 48 lattice at 3 = 6.2 discussed in the previous section. As mentioned above, 
the computations Avere performed using extended interpolating operators. The slope was 
determined for two values of the mass of the heavy quark (corresponding to Kh = 0.125, 0.133, 
slightly heavier and lighter, respectively, than the charm quark mass, [16]) and for a single 
value of the light quark mass (ki = 0.14144), which corresponds approximately to the mass 
of the strange quark [13]. 

(2) (2) 

We now discuss, in some detail, the evaluation of the ratios i?3 ' and R2 computed using 
extended interpolating operators. In particular, the discussion of the geometrical finite- 
volume effects of the previous section will have to be generalised to the case of extended 
interpolating operators^^. 

^^In order to obtain ^'(co) from the results for it is necessary to evaluate the derivative of /3{lo). This 
is a calculation in continuum perturbation theory and for itiq about the charm quark mass, ref. [8] found 
/?'(!) = -0.24. 

^^For R'i^\ even if local interpolating operators are used, eq. (22) can not be computed exactly unless the 
dependence of the matrix element M'^{n) on n is known. 



point (uj = 1) simply in terms of the ratios R'^'l In the limit L — )• 00, the derivative^^ of the 

(2) (2) 

Isgur-Wise function at a; = 1 is given by R^' — R2 , 




p(2) p(2) , p(2) 



(32) 
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(2) 

3.1 The ratio R2 (t) computed using Extended Interpolating Op- 
erators 

When extended interpolating operators are used, the wavefunction factor Z{p^) in eq. (8), 

depends on \p\ and therefore can not be talcen out of the sum over n. Thus, the exact 
12) 

computation of i?2 would require the knowledge of the wavefunction factors for all possible 
momenta of the lattice. However, in the limit L ^ 1, the dependence of on 
momentum can be approximated by: 



(33) 



(4) 



Z"{p') 



1 Z"{0) 

2 Z(0) ' 

d{p 



(34) 



(35) 



In addition, however, as was pointed out in the last section, discretisation errors could be 
significant for heavy-light mesons, and we saw how the changes suggested in eqs. (18-20) lead 
to a better understanding of the measured data. By including these corrections and using the 
approximation of eq. (33), i?2^^(t) for extended interpolating operators can be approximated 
by: 



R^2\t) 



at 



2 V4sinh^(^) sinh(ma) 



4i? 



(2) 



8cosh^(^) 



'2cosh2(^f ) ■ 



1 



Kt)sinh(^^ 



L 



c{t) sinh(^) 



(36) 



m 



hit) 

16sinh^(^) ' 4sinh^(^)sinh(ma) ' sinh^(ma) 



+ 



30^^ 



8R 



(2) 



4 sinh 



+ - 



at 



t^a^ 



sinh(ma) 



(37) 



(38) 



at 



4 sinh (^) sinh(ma) 



4i?(^) 



Z ) 



(39) 
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c(t) 



(40) 



60 - 40cosh2(^) 



X 




cosh2(^) sinh2(^) 



45 30 



+ 



a sinh (ma) 



180t 



(41) 



where we have neglected any contributions from excited states. 

The corresponding expression when the continuum energy-momentum dispersion relation is 
used, can be found by taking the limit a in eq. (36). 

Unlike the approximation used in deriving eq. (16), in obtaining eq. (36) the lattice spacing 
has been kept finite. We have performed exactly the sum over x^^ in eq. (8) and in calculating 
the sum over n, we have expanded the terms of the series in powers of yrn/L to second order 
(which requires an expansion of Z(n^p^j^) up to fourth order), and we have extended the 
limits of the sum between — oo to + oo. The accuracy of the different approximations 
involved in the derivation of eq. (36) can only be rigorously checked in the case of local 
interpolating operators. In that case, i?2 v) can be computed exactly and can be compared 
with the approximation of eq. (36) setting i?^^'*-^^ to zero. The results of this comparison 
are quite similar to those presented previously in Figs.l and 2, where possible finite-lattice- 
spacing effects were ignored^^. More specifically, for the physical situation which will be 
studied below (range of masses, ma = 0.6 — 0.85, range of times, t = 11 — 16 and for a 
spatial volume 24^) we have found that eq. (36) for local interpolating operators introduces 
an error of the order of 1.5% with respect to the exact result, computed by using eq. (14) 
and eqs.(18 - 20). 

In the case of extended interpolating operators, the validity of eq. (36) depends on the 
accuracy of the expansion of eq. (33), and Avill have to be checked by analysing the measured 
data. The numerical values of i?^'' and i?^^ will of course depend on the specific details 
of the smearing procedure. In this study we use gauge-invariant Jacobi smearing on the 
heavy-quark field, described in detail in ref. [19]. 

In principle, our preferred strategy is to constrain the mass of the meson to the value obtained 
from the behaviour of the heavy-light two-point function C2°\t) and to perform a two- 
parameter fit to the ratio R^"^ (t) in order to obtain i?^ ^ and i?^ ^ . However, for the expected 
range of values of i.e. i?^^ = 0(i?^^a^), i?2^^(t) is almost independent of i?^^ (-^z^ only 
appears in the volume correction term and its contribution is suppressed by powers of 1/t). 
Thus in practice, two-parameter fits prove to be very unstable; relatively small changes in 

^■*Eq.(36) for local interpolating operators and in the limit a ^ reduces to eq. (16). 
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Rfa-^ 




Kh = 0.125 


-2.75 


o 1 +1-4 
-1.7 


Kh = 0.133 


-2.38 


0.0 



Table 1: Results for i?^^'^^^ extracted from ordinary two-point functions fits. 

R^z^ can be compensated for by large ones in i?^^ . Because of the instability of the fits Avith 
two parameters, we modify our procedure as follows: we take as our central value i?^^ = 
and for i?^ ^^ke the value obtained from a single parameter fit to i?2 • To estimate 
the error in i?^^ we repeat the fits with other values of i?^^ as explained in the following 
paragraphs. A comparison of the predictions using this procedure with the measured data 
will serve as a further check of the accuracy of the several approximations performed. 

From the analysis of the usual heavy-light two-point functions {C'2\t)) we know the wave- 
function factors Z(0),Z(p^ = and Z{p^ = 2p^„-„) and therefore we can estimate the 
size of R^''''''^'^ ■ In Table 1 Ave present the values for the first and second deriA'atiA'es of the 
wavefunction factor Z obtained in this AA^ay. The results of Table 1, although affected by 
discretisation errors, show that in the conditions of the present study, the ratio R^^\r'^/R^p 
is indeed of order one for both heavy-quark masses considered, and one would expect a 
negligible dependence of R2'\t) on the exact value taken for H^^K 

In Table 2 we present the results for i?,^-* from different fits (including or neglecting correla- 
tions betAA-een different time shces) to R2'\t) using eq. (36), in all of which i?^^ has been fixed 
to zero^^. The meson masses have been constrained to the values obtained from measuring 
the usual heavy-light two-point functions. We see that the values obtained for i?^-* agree well 
with those quoted in Table 1. Moreover all the fits presented in Table 2 have small values 
of x^/dof, all of which supports the approximations used to account for the finite- volume 
corrections. In Figure 6, we shoAv the measured A^alues together Avith the fits using eq. (36) 
with R^z^ = 0, for both values of Kh. As can be seen in these plots, the agreement between 
data and our predictions is good and gives us confidence in our procedure. 

(2\ 

In the next subsection, we will use the values obtained for R^ here as input into our 
calculation of the slope of the Isgur-Wise function at zero-recoil. Fixing i?^^ to zero will 
introduce a small systematic error in our value of the slope, which will have to be estimated. 
Our strategy will be the following: i) For each k/j, the central value of the slope will be 
calculated by using the values of R^ obtained with the method outhned in this subsection 



^^For simplicity in the discussion, we do not present results obtained by using the continuum energy- 
momentum dispersion relation. In general, results obtained in this way reveal pronounced discrepancies 
between R^\t) and c'^\t), and have considerably higher values of x^/dof than those quoted in Table 2. 
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Kh = 0.125, R^z^a-^ = 


K,, = 0.133, R^z^a-^ = 




R^^^a-' 


ma 




ma 


corr. 

xVdof 


- 2.90 +^ 

— 5 

1.8/5 


[0.822 +3 


-2.51 +J 
5.5/5 


[0.598 


uncorr. 
xVdof 


- 2.91 +f 

— 5 

0.3/5 


[0.822 +5 


- 2.47 +^ 

— 6 

2.8/5 


[0.598 



Table 2: Results for i?^ for different fitting procedures. The fits were done for timeslices 
11, ... , 16. We present results obtained by using both correlated and uncorrelated fits. The 
masses have been obtained from C2'\t), and are shown in square brackets. 

(i.e. = 0). ii) In each case, we will recompute i?^^ and the slope of the Isgur-Wise 

function with two new fixed values of i?^^ given by three standard deviations from the 
central values of i?^^ quoted in Table 1. The spread of values obtained in this way will 
enable us to estimate the systematic uncertainties in our result for the slope associated with 
the approximations discussed in this subsection. 




Figure 6: Ratios R2 (t) for Kh = 0.125 and Kh = 0.133, together with the (uncorrelated) fits 
with masses constrained to their known values from ordinary two-point correlators. 
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(2) 

3.2 Extended interpolating operators : The ratio (t). 

The exact computation of (t), eq. (22), would require the knowledge of the wavefunction 
factors Z(n^p^j^) as well as the matrix elements M^in^p"^-^) for all the possible momenta of 
the lattice. Following the ideas of the last subsection, we will expand Z{n^p'^^-^) in powers 
of Pmin (eq. (33)), and similarly we will use an expansion for M*{n^p'^-^) of the type: 

M\nVmJ -M^(0) (l + (|) Vi?£) + (|) ) , (42) 
= ll^' ^^^^ 

with i?y defined in eq. (29). 

With the expansions of eqs. (33) and (42) and using the same kind of approximations as in 
the case of E!i\ can be approximated by 



2 V4sinh^(^) sinh(ma) ' ^ ' ) ^ 6'(t)sinh(^; 



a2 a2 /2cosh2(^f ) - 3 



- ('45') 

8cosh2(^) 24 V cosh2(^) y c'(t)sinh(^)' ^ ' 

where b'{t) and c'(t) are given by the expressions for h{t) and c{t) of eqs. (37-41) with the 
following substitutions: 

{{Rff + 2Rf) ^ {RfR^S + {Rf+R^^)l (46) 



2Rf ^ Rf+R^S- (47) 

In order to compute eq. (45), we set i?^'' to zero and we use the values for i?^' extracted 
from the fits to the ratio i',*), as in^ut. Thus, we have two free parametel, l^S and 
R^M ^ (involving the first and the second derivatives of the Isgur-Wise function at the zero- 
recoil point) which, in principle, should be determined from the fit to the ratio R\ (t). RjJ 
contains ^'(1), eq. (29), whereas R^^ is given by: 

R^S = ^(^'(1)--)- (48) 



^^In eq. (45) the "lattice" energy-momentum dispersion relation of eqs. (18- 20) has been used. The 
expression for the ordinary continuum energy-momentum dispersion relation could be found by taking the 
limit a in eq. (45). 
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The derivatives of the radiative corrections, which determine the relation between the deriva- 
tives of ^ and f at a; = 1, were calculated in [7] and [8], and around the mass of the charm 
quark (our k/^ = 0.125 = 0.133) quark is slightly heavier (lighter) than the charm quark, 
[16]) are approximately equal to = —0.24 and = 0.17. 

As in the case of Rf\t) and 7?.^-', for those values of ijj^-* and R^^j for which eq. (42) is a valid 
expansion, Rf\t) hardly depends on i?^-* and in the limit L — ^ oo is totally independent 
of it. On the other hand, we can get an estimate of the size of the ratio (-R^ ^ "^^/-^m ) 
by assuming a theoretical dependence of the Isgur-Wise function on a;, in a region close to 
Lo = 1. We have studied two different dependencies, widely used in the literature [20]: 



In both cases, for values of ~ 1.5 [5], we find that the ratio {R\J x m?/R)J) takes 
values of the order of 0.5. Therefore, the contribution of the R^^ term to the expansion of 
eq. (42) is suppressed by at least an extra factor with respect to the R^^ contribution, 
and furthermore the R^^ contribution to R^^\t) will be also suppressed by powers of 1/t. 
Because of the instability of the fits with two parameters we have followed two procedures: 

i) By using one of the u dependencies of eqs. (49-50), we can rewrite R^^^ in terms of 
RjJ and therefore perform just a one-parameter fit to the ratio B\ (t) in order to 
obtain R)J . We have checked that the difference in the final result for the slope of 
the Isgur-Wise function at ui = 1 depends only at the level of 1% on the particular ui 
dependence used. 

ii) Another possibility is to fix R^^ to zero; in this way i?^^ can again be extracted from a 
one-parameter fit to R^ (t). The results obtained in this way differ from the previous 
ones only at the level of 1-1.5% for Kh = 0.125 and at the level of 3% for Kh = 0.133 . 

To obtain our final results, we have decided to use method (ii). In this way, our estimate 
for the slope of the Isgur-Wise function will not depend on any particular choice of the 
theoretical dependence of the Isgur-Wise function on ui. The systematic error associated with 
this approximation, as we have discussed, could be of the order of 1.5% for Kh = 0.125 and 
3% for Kh = 0.133 , which is smaller than other sources of errors present in our calculation. 

In Table 3, we present our results for ^'{l). We quote numbers for ^'(1) obtained by fitting 

^2) (2) (2) 

Rl (t) to the functional form of eq. (45) and by fitting the difference of ratios i?3 (t) — -R2 (^) 
to the difference between the functions of eq. (45) and eq. (36). In all of the cases, i?^-* and 

(2') (2'] 

m have been obtained from our previous study of R2 (t) and C2 (t), respectively. The best 




(49) 




(50) 
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Kh = 0.125 


4^^ 


p(2) p(2) 


Kh = 0.133 




p(2) p(2) 
















corr. 

xVdof 


-2.90 +^ 

— 5 

1.8/5 


- 2.00 

— 15 

5.0/2 


-191 
0.7/2 


-2.51 +^ 
5.5/5 


15.9/2 


-1.61 
0.2/2 


uncorr. 
xVdof 


-2.91 

— 5 

0.3/5 


- 1 92 ^'^^ 

-L-^^ -14 

0.0/2 


- 1 94 '^^^ 

^■^^ -15 

0.1/2 


-2.47 

— 6 

2.8/5 


- 1.63 +^ 
0.4/2 


- 1.59 
0.0/2 



Table 3: Results for ^'(1) for different fitting procedures. The fits to the ratio of two-point 
correlators have been extracted from Table 2. The ratio of three-point correlators was fitted 
for timesUces 11, 12, 13, (in our simulation ty in eq. (22) is equal to 24). 



X^/dof are obtained when we look at the difference of ratios \t) —R2 (t)- Presumably, in 
this difference of ratios, the finite-volume corrections partially cancel and, therefore, results 
are less sensitive to the exact details of the procedure. 

As can be seen in Table 3 the results are rather stable under variations of the method used. 
For Kh = 0.125 a reasonable value for ^'(1) could be —1.92 ^j^. This is the vahie obtained 
using an uncorrelated fit with the mass of the meson constrained to its knoAvn value from 
ordinary two-point function fits. Correlated fits in principle give a more realistic estimate of 
the quality of the fit as measured by x^; however they are more sensitive to the uncertainties 
in the fitting function. The errors quoted for ^'(1) are only statistical and to them must be 
added those associated Avith our treatment of i?^'' and . As Avas discussed at the end of 
the last subsection, in order to estimate the uncertainties due to our ignorance of i?^"*- we 
have recomputed C'(l) foi" two more fixed values of R^^^ {R^^^ = — 3 and R^^^ = 6) obtained 
by allowing three standard deviations from the central value of i?^'' quoted in Table 1. The 
spread of values obtained in this way for ■^'(l) determines a further error of around of 1%- 
2%. This error, together with the 1.5% discussed before due to i?^ 5 gives a final estimate 
of our systematic uncertainties of the order of 3%. Thus, the value for ^'(1) quoted above, 
together with our estimate for the systematic uncertainties, leads to a result for the slope of 
the Isgur-Wise function at the zero-recoil point, for Kh = 0.125, of 

e'(l)=e'(l)-/5'(l) = -1.7 (51) 
where we have taken j3'{l) = -0.24 ([7]). 

(2) 

For Kh = 0.133, there is a greater uncertainty. In this case a correlated fit to i?3 has a 
large x^/dof, which may be an indication of a poorly-known fitting function. However, a 
correlated fit to the difference between ratios, R^^\t) — i?2^^(t), gives a good x^/dof and 
an answer for i^'(l) in excellent agreement with the answer from uncorrelated fits. Thus, 
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Figure 7: Ratios R]^ '{t) for = 0.125 and Kh = 0.133, together with the (uncorrelated) fits 
with the masses constrained to their known values from ordinary two-point correlators. 

following the same criterion as for Kh = 0.125, we quote a value for ^'(1) = (—1.63) or, 
equivalently, a value for the slope of the Isgur-Wise function at the zero-recoil point of 

e'(l)=e'(l)-/5'(l) = -1.4 (52) 

where we have included in a similar way to that above, the systematic errors associated with 
our choice of R^^ and R^m ^ ^'^d, for this value of the quark mass, we have taken = —0.20 
([7]). 

The values of eqs. (51-52) may be slightly higher than, but certainly compatible within errors 
with, the results of [7] given in the introduction, obtained with the traditional method of 
extrapolating from lu > 1 to lu = 1. 

(2,) 

In Figure 7 we show the ratio R^ '{t) together with the fit, for both values of Kh. 

By construction, the Isgur-Wise function, and therefore its derivatives, are independent of 
the mass of the heavy quark considered. Although within errors our results for the two heavy 
quarks studied here agree, the central values (eqs. (51-52)) may suggest some dependence 
on the heavy-quark mass. Trying to conclude that this effect is not merely statistical and 
trying to quantify the size of the 0{l/m,Q) corrections (derivatives at the zero-recoil point of 
7(0;), defined in eq. (24)) Avould require an extensive study, involving more than tAvo heavy 
quarks, which is outside the exploratory aim of this paper. Such a study is made, using the 
traditional method of extrapolating from a; > 1 to a; = 1, in ref. [7]. 
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Recently the CLEO collaboration has reported a value for ^'{l) = 1.01 ±0.15 ±0.09 obtained 
by looking at the semileptonic decay B D*, [21]. This number should not be compared 
directly Avith our results (eqs. (51-52)), as our results were obtained using a light antiquark 
heavier than the strange quark. However sum-rules ([22]) and quark-model ([23], [24]) calcu- 
lations suggest that the absolute value of the slope of the Isgur-Wise function at zero-recoil 
decreases with the mass of the light antiquark. 

4 Conclusions 

In this paper we have presented the results from an exploratory study implementing the 
method suggested in ref.[8] for computing the slope of the Isgur-Wise function at the zero- 
recoil point in lattice simulations. We have found large finite-volume corrections on currently 
available lattices, effects which will be present in all direct computations of derivatives of 
matrix elements with respect to external momenta. We have studied these corrections and 
have given their volume, mass and time dependence. We have seen that these finite-volume 
corrections do not have a dynamical origin, in contrast to those studied in [10]-[12], and the 
usual two- and three-point correlation functions (without the insertion) are not affected by 
them. Eq.(17) constitutes a simple way of estimating the expected size of these corrections 
for a given simulation. Furthermore, from this study, we have learnt how to account, in an 
approximate but sufficiently precise way, for these corrections and have been able to obtain 
reasonable values (^'(1) = —1-7 ^'{1) = —1.4 ^^) for the slope of the Isgur-Wise 

function, for two mesons composed of a heavy quark slightly heavier (ac/, = 0.125) and lighter 
{Kh = 0.133), respectively, than the charm quark, and in both light antiquark whose 

mass is about that of the strange quark (ki = 0.14144). These values compare well with 
those obtained in [7]. It would be interesting to study the dependence of the results on both 
heavy- and light-quark masses. 

We have seen that the understanding of the finite-volume effects is considerably simpler; 
indeed they seem to be fully under control, if local interpolating operators are used. The 
price for using local operators is that their overlap with the ground state is generally poorer, 
and one has to fit correlation functions at larger times in order to be able to neglect the 
contributions from excited states. Nevertheless, the calculation described in this paper 
should be repeated with local interpolating operators. 
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